Nanotube Quantum Dot Transport With Spin-Orbit Coupling and Interacting Leads 
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To date the impact of spin bias on the Kondo effect in nanotube quantum dots (QDs) is scarcely 
explored, especially in the case of significant Coulomb interaction between quantum dots and electric 
contact leads. Recent rapid experimental progress in nanotechnology opens new possibilities to 
study spin-bias-induced transport phenomena. Thus, of great interest are theoretical investigations 
of these transport properties in comparison with the case of a conventionally applied voltage for 
nanotube QDs interacting with contact leads. Such an investigation was carried out in this work 
where we analyzed the effects of a spin voltage as well as a conventionally applied voltage in a 
QD system with a different number of quantum states in the dot region in presence of Coulombic 
interaction between the quantum dot and two leads. The transport is described within the Keldysh 
non-equilibrium Green's function (NEGF) framework. We extended the NEGF treatment developed 
for noninteracting leads onto the case of four quantum states m = {a, A} = {±, ±} ( two spins and 
two energy subbands, which is the case for a nanotube QD) interacting with leads. Our derivation 
is based on the equation-of-motion technique and Langreth's theorem. For a Coulombic repulsion 
between the contacts and QD we obtain an expression for the current through QD for the four 
quantum states. To determine the parameters of the model Hamiltonian we used our previous 
calculations [1[ of the electronic properties of a symmetrical nanotube QD (5,5)/(10,0)_l/(5,5) in 
a tight binding model, where _1 denotes the length of the middle QD segment of a (10,0) zigzag 
nanotube. We calculated the density of electronic states with spin up and down for the case of a 
single QD without pseudospin states for an infinite Coulomb repulsion, in good agreement with the 
calculations of Yuan Li, et al. @]. Our calculation showed that the position of the conductance peaks 
nearest to zero is not affected by the strength of the QD-lead Coulombic interaction parameters. We 
also demonstrated that this interaction shifts the density of states to higher energies. The interplay 
between the Kondo effect and the bias is highly temperature-dependent and becomes significant 
only at low temperatures. Lastly, we found that the existence of four quantum states m — {a, A} = 
{±, ±} leads to abrupt changes in the density of states. In this case the values of the current 
are approximately ten times lower than for QD with only two quantum states m = {<r} = {±}. 
However, in the case of a conventional bias the current amplitudes in both cases are approximately 
the same. 



PACS numbers: 73.63.-b, 72.25.-b, 72.15.' 
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I. INTRODUCTION 



A fascinating many-body phenomenon in condensed 
matter is Kondo effect. It was extensively studied dur- 
ing the past years in connection to the increase of ap- 
plications of carbon nanotubes as quantum dot (QD) in 
variety nanoelectronic devices. Fabrications of carbon 
nanotubes opens new stage in the investigation of Kondo 
effect Q. In carbon nanotubes there is additional or- 
bital degrees of freedom which originates from two elec- 
tronic subbands near Fermi energy and plays a role of 
pseudospin. Recently, the results by Kuemmeth et al. [H 
demonstrate that the spin and orbital motion of electrons 
are strongly coupled, thereby breaking the SU(4) symme- 
try of electronic states in a such system. The spin-orbit 
coupling determines the filling order in the two-electron 
ground state [ij and leads to the fact that Kondo effects 
manifests as split resonant peaks in the differential con- 
ductance even at zero magnetic field. In recent years, ex- 
perimental progress opens new possibility for the study of 
spin-bias-induced transport in mesoscopic systems. For 
example spin bias can be achieved by controlling spin ac- 
cumulation at the biased contacts between ferromagnetic 



and nonmagnetic leads. This allows one to generate a 
pure spin current without any accompanying charge cur- 
rent [a, Q. However the effect of the spin bias on the 
Kondo effect remains weakly explored, especially in sit- 
uation where the Coulomb interaction of QD with leads 
are significant. More importantly, the Kondo effect in 
nanotube QD in the presence of Coulomb interactions 
between the QD and two adjacent leads has not been 
studied. The need of extending steady state transport 
formalism is required because new experimental results 
can't be explained with steady state transport formalism. 
Our goal was to derive such a formalism that could fit CN 
QD and can account for dot-lead electron interraction. 



II. MODEL AND FORMULATION 

We investigate model system with four-state QD with 
interplay between the spin and orbital (pseudospin) de- 
grees of freedom. This situation highly relevant to CN 
QD and is more general than two-state QD investigated 
previously (see for example [2|). We use the same strat- 
egy as in work [2| but the index that enumerates quan- 
tum dot state can run over four possible quantities that 
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lead sums and some another terms in formulas. There- 
fore this goal requires complete rederivation of formulas 
of . Of course the formulas of should be acquired 
as partial case from a more general expression. 

Let's use Keldysh formalism available for the case of 
one quantum dot with two energy levels split by spin- 
orbit coupling of electrons with different spin states @ 
and let's make it acceptable to the case of symmetrically 
connected carbon nanotubes (CNT) where there is addi- 
tional orbital degrees of freedom which originates from 
two electronic subbands near Fermi energy. Such a de- 
gree of freedom plays the role of pseudospin. For the 
nanotube quantum dot we adopted the following Hamil- 
tonian H = H c i + H c + Ht + H c d- The first term is the 
Hamiltonian which models an isolated CNT QD 
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where n m = d m d m and d m creates, d m annihilates elec- 
tron in the dot with m = {er, A} configuration, where 
{a, A} = {±, ±} are the spin and orbital quantum num- 
bers, U is the on-site Coulomb repulsion, e m = Ed — 
ctAAso/2 is the single-particle energy, its first term Sd 
is the basic dot level and the second accounts for spin- 
orbit coupling with the constant of a such interaction 
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Hamiltonian of 

the contacts, a — L,R - enumerates the left lead (L) 
and right (R) lead respectively, c\ km and c a k m cre- 
ates and annihilates electron with wave vector k and 
spin-pseudospin state m in lead a = L, R. The tun- 
neling between the dot and the leads is described by 
Ht = J2akm tam,Ca km dm + H.c. with spin and orbital 
conservation since the length of middle nanotube is very 
short (one translational period). In the case of conven- 
tional nanotube QD which is composed from one single 
walled carbon nanotube deposited on substrate to which 
gate voltage is applyed the role of the leads plays the 
same nanotube thus leads have the same orbital symme- 
try 0, H[ . Therefore for conventional NT QD spin and 
orbital conservation are valid. Our approach works for 
both cases conventional QD and QD from heterojunc- 
tions (n, n)/(2n, 0)_l/(n, n) but for obtaining numerical 
results we choose parameters for QD from heterojunc- 
tions (n, n)/(2n, 0)_l/(n, n). We also incorporate the 
Coulombic interaction between the charges on the quan- 
tum dot and two contacts. This interaction could play 
significant role especially due to the proximity of the con- 
tacts to dot. This last term in the Hamiltonian is given 
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Firstly we obtained a general formula of the current 
through CNT QD by considering the interaction between 
the QD and contacts. The current from the lead a 
through the barrier to the dot can be calculated from 



the time evolution of the occupation number operator of 
the lead a: J am = -e(N am ) = -^([H,N am ]) where 
N am = J2k c tkm c akm ■ Evaluating commutators we get 
ie ( \ 

Jam ~J~ f tam (Cakmdm) t am {d m Cakm) J ■ (3) 

By using the definition of less Green's function G< fcm = 
i( c akm(t')dm(t)) , the current can be written as 

Jam — — i?e[ 'y^ t am G akm (t, t)} . (4) 

k 

Let's consider the equation of motion for the time- 
ordered Green's function G^ fem (i,i') 

—ift-Q^G akm (t, t ) = e a k m G akm (t, t ) + t am G m (t, t ) — 



i/ afc \ T{d m (t)d+,(t')d m ,(t')c+ km (t')} 



(5) 

where we use the definition of time-ordered Greeen's 
function of dot G m (t,t') = -i(T{d m (t)d m {t')}), here T 
is the time ordering operator. The last term prevents 
obtaining a closed form for G akm (t,t'). At this step we 
need to use some approximations. The simplest possible 
solution is to use Hartree-Fock approximation 

-i{T{d m (t)d+, (t')dm> (t'K km (t')}) « 

n m ,(-i)(T{d m (t)c+ km (t')}) = nm>G a . kjn (t,t') 

where n m < = (d m ,d m >)- Therefore in Hartree-Fock ap- 
proximation the equation of motion for the time-ordered 
Green's function is: {-ih-§p - 4 fem )G^ fcm (t, t') = 
t^G^t,*'), where e' akm = e ak m + A ak , A a k = nl a k, 
n = J2 m n ™- If we d en °te Safcm(M') Green's function 
for operator {—ih-^j — e' a km) ano - operating them to pre- 
vious equation from right we get G^ fcm (i,i') = 
/ dtit am G t m {t,t 1 )g t akm (t 1 ,t'). Similarly the contour- 
ordered Green's function can be written as G a km(T, t') = 

J ' dTit* am G m (T,Tx)g akm [Ti,T'). Now using Langreth's 
theorem we can get G< km (t, t 1 ) 

GtkmM= I dr 1 t* am G m (T,T 1 )g< km (r 1 ,r')+ 



dritamGm(r,ri)g a akm(ri,r') (6) 



where the Green's function of leads for uncoupled system 
with modified energy levels e' akm is: #< fcm (*,*') = 
^nie^Je-^^'), g a a km {t,t') = i6(t> - 

t) e -< km (t-t') and f am {e' akm ) = (e '^'"" + l)" 1 is 
the Fermi distribution. After substitution of previous 
formulas to (6) and then to (4) we get 

%€. f 

Jam = 7TT / dzT am(z)(G m (z + A a (z)) + 

Znn J 

f am (z + A a {z)G m {z + A a (z))- 

fam(z + A a (z)G m (z + A a (z))) , 
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where T am (z) = 2nJ2k\ t am\ 2 S(z - e akm ) is the cou- 
pling function. We will consider proportional couplings 
to the leads Ti, m — x^Rm and symmetrize this formula 
by calculating the current J = xJl — (1 — x)Jr where 
J a = Y] m Jam and where the parameter x can be chosen 
as x = 1/(1 + x). For the case of proportional couplings 
we get: 

J = | X] / dzT m (z)(f Lm (z + A L (z))p m {z + A L {z))- 

m 

fRm(z + A R (z))p m (z + A R (z)) , 

(7) 

where T m {z) = ^S^lty P™{z) = -±Im[G r m (z)]. 
Using last formula only demands to know G r m (z). With 
the goal of finding retarded Green's function G r m (z) we 
use standard equation of motion technique based on the 
general relation 

< F 1 \F 2 > r z = {{Fi\F 2 })+ < [Ft,H]\F 2 » r . 

Writing this relation for the case Fi — d m and F 2 = e2+ 
we get 

?#P 

< d p \d+ > (z - e p ) = 1 + < d+d q d p \d+ > + 

i 

X]*ap < c a k P \d+ > -2^t* ap < c akp d+d p \d+ > + 

ak ak 

Y, Iak ^ c ikq c *k q d P \d+ > . 

akq 

(8) 

In this formula and below under by <C F\ , F 2 ^> we un- 
derstand <C F\,F 2 ^ r for the sake of simpli?cation in 
writing formulas. For the dot region we use indexes p, q, r 
and a, j3 for numbering the left a — L, j3 — L and right 
a = R, /3 = R contacts. For wave vector in leads we 
use k,k . In the right hand side of expression (8) appears 
new Green's function for each of which we also can write 
equations of motion 

< c Qfcp |d+ > (z - e akp ) = t ap < d p \d+ > - 

> . 

In the last formula we neglect three-operator terms be- 
cause they contain high order spin correlation. We do the 
same for all other formulas. Whenever a term with F\ 
contains five or more operators in -C F\\d p 3> we drop 
it. For p 7^ q we can write: 

<C d+d q d p \d p » (z — Ep — U) = n q + 

C P « c akp d+d q \d+ » + ^ « c akq d+d p \d+ > - 

ak ak 

J2 ta 1 < C tkq d q d p\ d t » ' 
ak 

(9) 



< c akp d q ~d q \d p v » (z - e akp - I ak ) = 
i Qp < d+d q d p \d+ > +y^tg g < c akp C[3 k > q d+\d+ > - 

/3fc' 
/3/c' 

(10) 

^ Cakpd p d p \d p (z ^akp Iak) — 
^ ' ^/3p Cf3k'pC a kpd^\d p 3> — 

^ t/9 P < c+ klp c akp d p \d+ > , (11) 

/3/c' 

< c akq d+d p \d+ » (z - e p + e q - e akp - I a k) = (cakqd+) + 

t aq < d+d q d p \d+ > + < CctkqCfik' 

k' 

< C pk'q C »kqdp\d+ » ' 

/3fc' 

(12) 

< 4fc 9 d 9 d pMp > (z - C/ - £p - e 9 + e Qfcp ) = 

XI l }l < C akq C l3k'qd p \d+ » + 
/3/c' 

< c tk q c l3k'pd q \d+ » -t* g < d+d q d p \d+ > , 

/3/c' 

(13) 

^ c ife 9 d ^pMp > - £ P + £ g + e Q fe 9 + ^afc) = 

/3/c' 

X */s« « 4*'«4fc,^K > • ( 14 ) 

/3/c' 

We used decoupling approximation based on the follow- 
ing rules [13 (YX) = 0, < FX 1 X 2 |d+ >« (XxXa) 

Y\d p ^> where X represents an operator of contacts 
and Y an operator of quantum dot. 

<C dp|d+ > (Z - £ p - X I a k{cl kq C a kq)) = 
akq 

q¥=p 

1 + < d+d q d p \d+ > + 

X*«P ^ C "fepMp > ~ 2 X^"P ^ C akpd p d p \d+ > , 
q/c a/c 

(15) 

< c a kp\d p » (z - e Qfc p) = t Qp (l - 2n Qfe p) <C » 

«7^P 

~^~Iakp Cakpd p d p \d p A-I a k ^ ^ 

d+d q \d+ > , 
(16) 
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< c a k P dpd p \d+ > (z - e akp - I ak ) = 

«d p |rf+>, (17) 

<C d q d q d p \d p » (z — e p — U) = n q + 

Y,K P ^c akp d+d q \d+^> + 

ak 

t* aq « C akq d+d p \d+ » -t^ < C+ kq d q d p \d+ ») , (18) 

< Cafepd+dqldp > (z - e Q/i; p - I Q fc) = 

tap < d, d q d p \d+ > , (19) 

< c a kpd q d p \d+ > (z - e p + e q - e akq - I ak ) = 

t aq <C d+d q d p \d+ > -t aq n akq < > , (20) 

^ c tk q d i d p\ d t y> (z - U - e p - e q + e akq ) = 
t* aq n akq <C > -t* aq < d+d q d p \d+ > . (21) 

From (16) using (18, 19) we get: 

« c akp \d+ »= tap(1 ~ 2 " afcp) « d„|< » - 



z-e. 



akp 



lak^c 



{z £akp)(,Z ^akp Iak) 

q¥=p 

(n akp < d p \d+ > < d+d q d p \d+ >) . 

The other formulas are obtained in a similar way: 
< dp|d+ > f ^ - £ P - ^I a k{cl kq c akq ) 

£ 



otkq 



\tap\ ^1 + lak^akp -j j 



^akp Z tctkp Iak 



a fe 



(•Z £akp){z ^akp Iak) 

q=£p 

£ < d+d q d p \d+ > 



After introduction 



It 12 



4 0) (*) = £r3? 



^ ^akp 



4^) = £7- 



a/c 



(•2 e a£;p)(^ ^akp Iak) 



l^-ap | lak^akp 



^ ^afcp) ^akp Iak ) 

E P 2) = £ Walt, - ( 22 ) 



we get: 

« d p < » (z-Sp- E< 2 > - Ef(z) - Ef( Z )) 



+ • (23) 



\t 2 

^ap 



From (18) using (19,20,21) in a similar way we get: 
< d+d q d p \d+ > ( z - £ p - t/ - £ 

-£( 



^ ^akp -^-ak 



L aq | 



Z £p -\~ £q ^akq Iak 



+ 



L aq | 



Z U Sp Sq + ^akq 



i g -«dpM+»£(— 



\taq\ Tic 



ak 



Z Sp ~\~ S q ^akq Iak 

\^aq\ Tlakq 1 
Z U Sp Sq ~\~ ^akq 



+ 



After introducing: 

^(*) = E( |/n " 12 



Z £p -f- Sq ^akq Iak 

l^agj ^ 

Z £"p Sq -\~ Cctkq 



-+ 



^4 ^ _ |^ag|_^afcg ^ 

takq Iak 



ak 



\tgq\ n akq \ 



Z - U - S p - S q + 

^akq 

we can rewrite it as: 

<C d+d q d p \d+ » (z-Sp-U - £ 



|^ap I 



^ ^akp Iak 



A q p(z)) = n q - < > i 9 p(z) . 

(25) 

The last formula (25) allows as to write (23) in the form: 
« d p \d+ » (z-e p - Sf^) - SW^) - S< 2 ) + 

(u+^(z))y: 

q 

q¥^i 



, z-ep-U -Y>f\z)-A qp {z) 
q¥=p 



q z Ep U Sp ^ i z ) A q p(z) 



(26) 



where t p 0) (z) = J2 a k z-ll7-i ak ■ The last e q uation (26) 
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can be solved for <C d p \dp ^>= G p (z) : 



G p (z) = 1 



,, z — Sp — U — '(z) — v4,jp(z) 
/( Z -s r Sf(z)^W(z)-S( 2 ) + 



(27) 



here Tlakp — fapi^akp^i fap^akp) (cxp( , r 

l) 1 - Fermi distribution, - chemical potential of 
lead a for spin and orbital quantum number p = {a, A}. 

Let us explore the limit U — > oo. By taking this limit 
in (27) we get 



G p (z) 



i-E 



q U <1 



z-e p - Z p 0) (z) - t p L> (z) - E( 2 ) - S p Jj (z) 



-.(3), 



Green's function G p (z) depends on unknown quantities 
n g = (dqd q ). But n g = — ^ J dzG^(z) and using Lan- 
greth's theorem we get 

w d /t ,(^(z) + r^(z)/H,(z) fm[ j 
ttj r L? (z) + r fl? (z) 9 

(29) 

Formulas (27) (or in the case U —y oo approximation 
(28)) with (29) allows us to self-consistently compute n q 
and G p (z). The last thing needed for using this pro- 
cedure is to rewrite formulas (22,24) in integral form 
using the definition of the coupling function T am (z) = 
2lr J2k l'amP*( z - £afcm)- For example 



\tgq\ n akq 



(28) 
The 



2- takp 



2tt 



(30) 



Therefore by choosing a reasonable predefined form of 
coupling function r am (z) we can self-consistency solve 
equation (28,29) and find the retarded Green function 
G p {z + irj) by using complex energy argument with suf- 
ficiently small imaginary part 77 — > + and further use it 
for calculating the current by formula (7). 



III. NUMERICAL RESULTS 

We now apply the obtained formulas for nonequilib- 
rium transport in nanotube quantum dot with the cen- 
tral part coupled to interacting nanotube leads. In our 



numerical calculations we assumed that nanotube quan- 
tum dot is symmetrically coupled to contacts from nan- 
otubes with Lorentzian linewidth 2D, Tl p (z) — Tr p (z) = 

jjrrp ■ From formula (30) it is obvious that all energy 
scale is dictated by the coupling function 7. In other 
words setting 7 = 1 means that the energy is measured in 
7 units. For comparison with experiments and other the- 
oretical works the resulted quantities can be transformed 
in electron Volt energy scale by adopting 7 = lOmeV [ill ]. 
For example the energy cutoff parameters D = 500 
used in our computation means D = 5OO7 = 5eV. 
The dimensionless temperature T = 2 corresponds to 
T = 27^7 = 20mA'. The contact dot interaction we 
used a flat-band profile I a (z) — i a 8(D — \z\), where 
%L,iR are parameters that describe Coulomb interaction 
with left and right leads respectively. We next chose 
parameters which concerns the nanotube dot region. 
For nanotube quantum dot we adopted the dependence 
£m = £ {a,\} — £ d — cAAso/2, here e c i is the basic dot 
level. The value of e<j was taken from our previous calcu- 
lation [l[ which accounts for the geometry of the NTQD 
(5, 5)/(10, 0)J/(5, 5) and depends on the length I (here 
I = 1 one unit cell) of the central fragment. By varying 
the length, different set of peaks appear in the density of 
states of the fragment. The tight binding approach used 
in this work did not include accounts electron-electron 
interactions. The present work does so, but it does 
not include the actual atomic structure of the nanotube 
system. We only took one peak closest to the Fermi 
level and set Ed to its energy. For (5, 5)/(10, 0)_l/(5, 5) 
nanotube quantum dot Ed — 0.006£c-c was calculated 
by us n| . Using for graphene hopping parameter the 
value [l2j tc-c ~ 2.66eV we get in dimensionless units 
£d = lhmeV = I.57. The model used does not either 
consider spatial spin correlation which is taken into ac- 
count through some sort of mean field approximation 
through Age?, in which the main contribution gives 
curvature [13j coupling A^q' 11 m 1.6meV/d[nm]. Di- 
ameter of (10,0) nanotube is d = 3.9nm which gives 
A so = 0.2meV = 0.02 7 . 

Applying spin bias V s means ^l\x = ~^lix = V s /2 
and [Irx\ — ~/j.rix — —V s /2. By swap indexes L with R 
or which is the same assigning hl^x = —^lix = —V s /2, 
MiffA — ~^Rix — Vs/2- It is possible to mirror this 
dependence by OY axis but we decided to leave it un- 
changed for reason of distinguishing between graphs for 
applied spin and conventional bias on basis of slope of a 
such dependences. 

In this work we used U — > 00 approximation. Studying 
the effect of U will be devoted to another work. First (see 
Fig. 1) we calculated the current- voltage characteristics 
for spin bias at T = 1 temperature for different Coulomb 
coupling to contacts jj, = in = green, ii = ir = 0.05 
red curve. We also plotted its derivative (see Fig. 2), 
the so called differential conductance, as a function of 
applied spin bias V s . In all graphs dimensionless quan- 
tities are used. To get dimension quantities we need to 
take 7 value (for example 7 = lOmeU). The position of 
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voltage (y/e=10 mV) 

Fig. 1. Current for the case of spin bias voltage. Two states 
m = {a} QD at temperature T=l without QD-lead interaction 
i£ = in = (green line) and with «i, = ijj = 0.05 (red line). 

peaks on Fig 1. near 0.1 in dimensionless units means 
0.1 * 10^ = ImV. On the Current- Voltage charac- 
teristics there is a gap which is connected with Coulomb 
blockade regime. As we can see from Fig. 3 calculated for 
higher temperature T = 100, at this temperature the gap 
is already washed out. At low temperatures the differen- 
tial conductance (Fig. 2) exhibits narrow peaks and the 
nearest to zero voltage peaks are sufficiently extended in 
the case of presence of Coulomb interactions with leads 
ih = iR = 0.05. 

As we can see on Fig. 3,4 the temperature in- 
crease blurs the peaks of the differential conductance 
and smooth the current-voltage dependence. Therefore 
the interplay between the Kondo effect and the bias is 
highly temperature-dependent and is significant only at 
low temperature. 

We also analyzed the density of states of quantum 
dot p p (z) — —(X/ir)ImGp(z) at applied spin-bias voltage 
V s = 0.3. Fig. 5 depicts the calculated density of states 
for spin down electrons for different Coulomb QD-leads 
interaction. The absence of interaction (ii = ip_ = 0) is 
depicted by green curve, the presence (i^ = ir = 0.05) 
plotted by the red curve. In both cases voltage V s = 0.3 
is applied. As we can see there are two Kondo peaks 
which are split from each other due to the applied Volt- 
age. When V s = our calculation gives one Kondo peak. 
These results are in good agreement with 0]. 

As we can see from Fig. 5, the Coulomb dot-lead inter- 
action shifts the density of states in high energy region. 
As we can see the position of nearest to zero conductance 
peaks are not affected by the strength of the Coulomb 
dot-lead interaction. 

The next step was to study the nanotube QD in which 
four states m — {cr, A} = {±, ±} exist. In this case it is 
necessary to sum from q=l to 4 in formulas like (28). As 
earlier we consider spin bias voltage. Fig. 6 depicts calcu- 
lated density of states for four different quantum states 
m = {a, A} = {±, ±}. As we can see the existence of 
four quantum states leads to abrupt changes in the den- 
sity of states (see Fig. 5,6). In Fig. 7 the current-voltage 
characteristic for applied spin-bias voltage is presented. 
As we can see this dependences is about 10 times less 
than that for the case of two states m = {a} = {±} in 
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Fig. 2. Conductance for the case of spin bias voltage. Two states 
m = {<t} QD at temperature T = 1 without QD-lead interaction 
ijj = iji = (green) and with i^ = in = 0.05 (red line) 

quantum dot. The second main difference in differential 
conductance is the presence of a second pair of strong 
peaks in the case ii, = in = 0.05 which is absent in the 
case ih = ir : = 0. 

Lastly, we investigated the effect of conventional bias, 
namely hl^x = A«Z4A = -V/2 and n m \ = ^ Rix = V/2. 
The same parameters (Fig. 1,2) as in the case of conven- 
tional bias (Fig. 9,10) give about hundred times lager 
current and lead to a different structure of peaks of con- 
ductance. 

The conventional bias in the case of four states m — 
{a, A} = {±, ±} (nanotube quantum dot) gives (Fig. 
11,12) about three hundred times lager current and yields 
a different structure of peaks in the differential conduc- 
tance. As we can see (Fig. 8,10,12) the presence of QD- 
lead interaction leads to formation of new pair of peaks. 
The present results allows us to conclude that the QD- 
lead interaction is important and in all cases where it is 
strong it must be considered. 



IV. CONCLUSIONS 

We derived generalized Kubo formula for the com- 
putation of conductivity NTQD and formulas for self- 
consistent evaluation Green's function of NTQD in pres- 
ence of Coulomb repulsion of electrons in QD region and 
electrons in lead. We show that presence of four states 
in NTQD with m = {a, A} = {±,±} leads to abrupt 
changes in the density of electronic states and to lower- 
ing the current amplitude by approximately a factor of 
ten compared to QD with only two states m = {a} = {±} 
(without pseudospin). But in case of conventional bias, 
current amplitudes are approximately the same. We also 
show that the Coulomb interaction ii and ir with the left 
and right leads shifts the QD electron density of states to 
high energy region and the interplay between the Kondo 
effect and the bias is highly temperature-dependent and 
becomes significant only at low temperatures. Taking 
into account Coulomb dot-lead interaction leads to the 
emergence of a second pair of peaks in the differential 
conductance with significant amplitude. We also show 
that a such an interaction does not shift the pair of peaks 
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voltage (y/e) 

Fig. 3. Current for the case of spin bias voltage. Two states 
m = {cr} QD at temperature T = 100 without QD-lead interaction 
= in = (green) and with i^ = in = 0.05 (red) 
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Fig. 4- Conductance for the case of spin bias voltage. Two states 
m = {cr} at temperature T = 100 without QD-lead interaction 
ijj = in = (green) and with i^ = in = 0.05 (red) 





energy (y) 

Fig. 5. Density of states for QD with two quantum states 
m = {cr}, a spin voltage of v s = 0.3, spin down as well as without 
(green line) and with (red line) a QD-lead Coulombic interaction 



energy (y) 

Fig. 6. Density of states for QD with four quantum states 
i = {cr, A} and a spin voltage of V s = 0.3 when m = {+, +} or 
{ — , — } green) and m = {+, — } or m = { — , +} (red line). 
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Fig. 7. Current for the case of spin-bias voltage. Four states 
= {a, A} QD at temperature T=l, without QD-lead interaction 
i^ = in = (green) and with i^ = in = 0.05 (red line) 
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Fig. 8. Conductance for the case of spin-bias voltage. Four states 
m = {cr, A} QD at temperature T=l, without QD-lead interaction 
ilj = in = (green) and with i^ = in = 0.05 (red line) 




-0.6 -0.4 -0.2 0.2 0.4 0.6 



voltage (y/e) 

Fig. 9. Current for the case of conventional voltage. Two states 
m = {cr} QD at temperature T=l without QD-lead interaction 
iL = in — (green) and with i^ = in = 0.05 (red line) 
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Fig. 10. Conductance for the case of conventional voltage. Two 
states m = {cr} QD at temperature T=l without QD-lead 
interaction i^ = in = (green) and with i^ = in = 0.05 (red line) 
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voltage (y/e) 

Fig. 11. Current for the case of conventional voltage. Four states 
m = {a, A} QD at temperature T=l, without QD-lead interaction 
lL = ir = (green) and with %l = in = 0.05 (red line) 




voltage (y/e) 

Fig. 12. Conductance for the case of conventional voltage. Four 
states m = {<r, A} QD at temperature T=l, without QD-lead 
interaction «l = *B = (green) and with ih = in = 0.05 (red line) 



(first pair) nearest to the Fermi level. 
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